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Abstract 

Similarly to the standard effective range expansion that is done near the threshold 
energy, we obtain a generalized power-series expansion of the multi-channel Jost- 
matrix that can be done near an arbitrary point on the Riemann surface of the 
energy within the domain of its analyticity. In order to do this, we analytically 
factorize its momentum dependencies at all the branching points on the Riemann 
surface. The remaining single-valued matrix functions of the energy are then ex- 
panded in the power-series near an arbitrary point in the domain of the complex 
energy plane where it is analytic. A systematic and accurate procedure has been 
developed for calculating the expansion coefficients. This means that near an ar- 
bitrary point in the domain of physically interesting complex energies it is possible 
to obtain a semi-analytic expression for the Jost-matrix (and therefore for the S- 
matrix) and use it, for example, to locate the spectral points (bound and resonant 
states) as the S-matrix poles. 
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1 Introduction 

Taylor-type power-series expansions are very common in physics. In quantum scattering 
theory, the most frequently used expansion of this kind is known as the effective-range 
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expansion. In the case of a short-range potential, it represents the cotangent function of 
the scattering phase shift Sg(k) in the form 

oo 

k 2i+l COt 6 e (k) = Y, C tnk 2n ■ (1) 

n=0 

Here the right-hand side is a sum of terms proportional to even powers of the collision 
momentum k, and cg n are energy independent expansion coefficients. Originally (see, for 
example, Ref. [1] and Ref. [2] for a historical review), this expansion was suggested in 
nuclear physics for the S'-wave nucleon-nucleon scattering in the form 

k cot 5 (k) = -- + K k 2 - Pr 3 k 4 + Qr 5 k 6 + • • • , (2) 

where the first two parameters on the right-hand side, namely, a and ro, were called the 
scattering length and the effective radius. The parameters in the higher terms of this 
expansion (P, Q, etc.) were known as the shape parameters. This name originated from 
their dependence on the shape of a chosen iViV-potential, while a and ro were mainly 
determined by the strength of the interaction and the range of the force. 

Originally, some rather inconvenient integral formulae (see, for example, Ref. [2]) 
for calculating the parameters on the right hand side of Eq. (j2j) were suggested. The 
expansion, however, was mainly used to parametrize the experimental data by choosing 
appropriate values for the first few coefficients. It is still frequently used for such purposes 
not only in nuclear physics but also for parametrizing the low-energy collisions between 
atoms and molecules (see, for example, Refs. [31 HJ [7J [8]). The expansion is also 
useful for simplified analytical description of the scattering near the threshold energy. For 
example, in one of the recent works of this type the effective-range expansion is used to 
study the causality constraints on low-energy universality [9]. 

A workable and rather accurate method for calculating the low-energy parameters a, 
ro, P, etc., was developed about fifty years ago [TOj HU [121 1131 EE] within the variable- 
phase approach. Recently the same equations were re-derived in Ref. [T5] using the 
technique of canonical quantization. In this method, the low-energy parameters are the 
asymptotic values (r — > oo) of the corresponding radial functions a(r), r (r), P(r), etc., 
obeying a system of first-order differential equations. At a fixed r, these functions give 
the corresponding low-energy parameters for the potential which is cut off at this value 
of the radius. 



2 



The main drawback of using the differential equations for a(r), r (r), P(r), etc., is 
that such a method is inconvenient in the case of potentials supporting bound states. 
Indeed, with the radius increasing from zero to infinity, the potential grows from nothing 
to its complete shape and its bound states appear one by one at certain values of r. Each 
new bound state initially appears with zero binding energy which corresponds to infinite 
scattering length. Therefore the function a(r) has singularities, which cause serious prob- 
lems when the system of differential equations is solved numerically. This difficulty is 
avoided if instead of a(r), ro(r), P{r), etc., we consider the expansion coefficients of the 
Jost function for the potential cut off at r. These coefficients also obey a system of differ- 
ential equations but do not have any singularities. The standard low-energy parameters 
are obtained from them via simple algebra |16| . 

The traditional approach has one significant limitation, namely, the effective-range 
expansion is only applicable near the point k — 0, i.e. when the energy is close to the 
threshold. In Ref. [16], the expansion (0Q) is generalized by writing it as a series of pow- 
ers of (k — k ) with an arbitrary complex k . By doing this, it is possible, for example, 
to explore the complex fc-plane in search for resonances. When fc 7^ 0, the generalized 
effective-range expansion is no longer a low-energy approximation. In other words, in- 
stead of making the expansion around the threshold, one can do it around any point in 
the complex plane. 

Yet another limitation of the original expansion ([T]) is that it is formulated for a 
single- channel problem. For multi-channel problems, something similar was also suggested 
[T7J [HI [191 [20], |21j [22] . In these papers, the authors tried to keep as much similarity with 
the single-channel case as possible. Their adherence to the functions of the type cot 5 
limited the flexibility and clarity of the suggested equations. 

In the present paper, we use a different approach in generalizing the effective-range 
expansion to multi-channel problems. Instead of considering the channel phase shifts and 
their cotangents, we look at the problem from a more general point of view. As a first 
step, we analyze the structure of the Jost matrix and explicitly factorize the dependencies 
of its matrix elements on odd powers of the channel momenta. These factors determine 
all the branching points of the Riemann surface of the energy, while the remaining factors 
are single- valued holomorphic functions defined on a simple energy plane. These functions 
are then expanded in power series of (E — Eq), where Eq is an arbitrary complex energy. 
As a result, we obtain an expression for the Jost matrix such that each of its elements 
is a product of a nonanalytic "branching" factor and a power series. In this way, we are 
not limited to the threshold points, but expand the Jost matrix practically anywhere on 
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the Riemann surface. When the Jost matrix is obtained in such a semi-analytic form 
(the first several terms of the expansion), the S matrix and all the observables can be 
easily calculated. Although it is possible, there is no need in introducing the generalized 
scattering length and other parameters. The Jost matrix expansion coefficients are more 
simple, clear, and convenient for this purpose. 



2 Mult i- channel Schrodinger problem 

Consider a quantum mechanical two-body problem which, after separation of the motion 
of its center of mass, is reduced to an effective problem of one body whose dynamics is 
governed by the Hamiltonian 

H = H +U + h (3) 

consisting of the free-motion part H Q , the interaction operator U, and the Hamiltonian h 
that describes the internal dynamics in the moving body (for example, the internal states 
of colliding atoms). Each internal state of the body corresponds to a different channel of 
the scattering process. In general, there are infinite number of the internal states, i.e. the 
eigenstates of h, 

h\n) = E n \n), n = 1, 2, 3, . . . oo (4) 

We assume that only N internal states are important. In other words, we approximate 
the internal Hamiltonian by the N terms 

TV 

h ^2 \ n )En(n\ . (5) 

The total Hamiltonian taken in the representation of the relative coordinate f and sand- 
wiched between (n\ and \n') ) becomes the N x N matrix 

h 2 

H nn > = —S nn '- — A? + U nn >(r) + E n 5 nn > , (6) 

where /z n is the reduced mass in the channel n. Therefore the eigenvalue equation if \l/ = 
E^/ is reduced to a set of coupled differential equations for the channel wave functions 

N 

ME, r) = U nn ,{f)^{E, r) . (7) 

n'=l 



2^ 



■A r + (E — E n 
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The eigenenergies E n of the internal Hamiltonian h are the thresholds for the correspond- 
ing channels. An eigenstate of the Hamiltonian (jH]), corresponding to the eigenvalue E, is 
a column-matrix 

( ME,r) \ 



\ *p N {E,r) J 

where each row describes the relative motion in the corresponding channel. 



(8) 



To avoid unnecessary complicated notation, we assume that the relative motion in 
each channel has a single value i n of the angular momentum. This does not compromise 
the generality of our consideration. Indeed, we can always treat the states with different 
values of angular momentum (or of any other quantum number) as different channels, but 
perhaps with the same threshold energies. 



Therefore the angular dependence of the channel wave function can be factorized in 
the simple way 

11 f R rl 

(9) 



This reduces Eq. © to the system of coupled equations for the radial parts of the channel 
wave functions, 



A' 



(r)u n >(E,f) 



n'=l 



where 



V nn > (r) 



2/i n 

h 2 



Y LmJ&i ¥)Unn'(r)Y en , mn , (9, if) dQ f 



and the channel momenta are defined as 



2 fi r , 
ft? 



{E - E n ) . 



(10) 



111) 



(12) 



In what follows, we assume that the interaction potential is non-singular and of short 
range, i.e. that all its matrix elements are less singular than 1/r 2 at the origin and expo- 
nentially vanish at infinity. In principle, the class of acceptable potentials can be much 
wider. This however would require a complicated analysis of the analytic properties of the 
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multi-channel Jost matrix in order to identify the domain of complex energies where the 
Jost matrix can be analytically continued to. With the exponentially decaying potentials, 
we can rely on such an analysis given in Ref. [23J where the domain of analyticity of the 
multi-channel T-matrix and S-matrix is rigorously established. Therefore by narrowing 
the class of the potentials, we sacrifice the generality for the sake of clarity. 

The boundary conditions for Eqs. ffTUl) are derived from the requirement that any 
physical solution must be regular at the point r = and have special behaviour when 
r — > oo, which is different for bound, resonant, and scattering states. In Sect. IU we 
consider these conditions in more detail. 



3 Multi-channel Jost matrix 



A system of N linear second-order differential equations of the type (flOl) has 2N linearly 
independent column-solutions and only half of them are regular at the origin (see, for 
example, Ref. [21] )• Combining these regular columns in a square matrix, we obtain the 
so called fundamental matrix of the regular solutions, 



ii 2 (£,r) 
^22 (E,r) 



4> 2N (E,r] 



\(f> m (E,r) <j) N 2(E } r) 



>NN 



(13) 



(E,r)J 



Any other regular solution can only be a linear combination of the columns of this funda- 
mental matrix. In particular, any physical solution of Eqs. ( TTOT) is a linear combination 
of the columns of the corresponding fundamental matrix, 



\u N J 



Ci 



+ c 2 



>1\ 



\0NlJ 



'12 
i>22 



\ 



/ 4>in\ 



+ --- + C 



N 



\07V2/ 



>2N 



(14) 



\4>nn) 



This guaranties its correct behaviour when r — » 0. As far as the asymptotic behaviour 
(r oo) is concerned, its correct form can be achieved by proper choice of the combina- 
tion coefficients C n . 
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Since the potential is of a short range, far away from the origin the right hand sides 
of Eqs. flTUl) vanish and these equations decouple, 



d r -\- k 



u n (E,r) 







when r — > oo . 



(15) 



These are the Riccati-Bessel equations. As a pair of linearly independent solutions of each 
of them, we can take the Riccati-Hankel functions h^(k n r). 



Despite the fact that equations f|T5|) are decoupled, we can still treat them as a system. 
It therefore has 2N linearly independent column- solutions which can be chosen in many 
different ways. The most convenient choice is the following set of 2N columns grouped in 
two square matrices 



\far) 







jy(out) 



/ h 



\ 



(+) 



(far) 











h£\far) 









\ 



ht\k N r)J 



(16) 








h[t\k N r)/ 



(17) 



that represent the in-coming and out-going spherical waves in all the channels. These 
2N columns form a basis in the space of solutions. In other words, any column-solution 
of Eq. ( !T5|) is a linear combination of these 2N columns. In particular, each column of 
matrix ( Tl3l) at large distances becomes such a combination. The combination coefficients 
have two subscripts: one to indicate which column of (fl3l) is expanded and the other is 
the summation subscript. Similarly to Eqs. ffT6| [T7|) . we can group these coefficients in 
square matrices (they depend on the choice of the energy, but do not depend on r). Thus 
we have 

$(E,r) — > W {[n) (E,r)F {[n \E) + W (out) (E,r)F {ont) (E) , (18) 

1 — >oo 

where, by analogy with the single-channel case (see, for example, Ref. [22]), the energy- 
dependent matrices F^ n ^ out \E) can be called Jost matrices. It is not difficult to show 
that they determine the S'-matrix 



S(E) = F (out) (£) [F^(E 



(19) 
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and thus give complete description of the underlying physical system. The spectral points 
E = S n (bound states and resonances) are those where the inverse matrix [F^(E)] 
does not exist, i.e. the points where 

detF( [n \£ n ) = . (20) 



4 Transformation of the Schrodinger equation 



At large distances, the fundamental regular matrix $>(E,r) is a linear combination (|18p 
with the r-independent coefficient matrices F^ n ' ont ^(E). We can, however, look for <&(E, r) 
in the same form at any point r, but with the coefficient matrices depending on r, 



${E,r) = W (in) (E,r)F (in \E,r) + W {out) (E,r)^ out) (E,> 



(21) 



Now, instead of one unknown function <&(E, r), we have two unknown functions, ^^(E, r) 
and J-(° ut \E,r), which therefore cannot be independent of each other. In principle, we 
can arbitrarily impose any (reasonable) additional condition relating them. The most 
convenient is to demand that 



W (in) (E, r ) ^ J- (in) (E,r) + W (out) (E, r) ^ (out) (E,r) 
or or 



, 



(22) 



which is standard in the theory of differential equations, and is known as the Lagrange 
condition within the variation parameters method [26J. Considering Eq. (Tl8|) . we see that 
the Lagrange condition is certainly satisfied at large distances. Therefore, by imposing it, 
we do not change the asymptotic behaviour of the solution. 



Starting from the coupled-channel radial Schrodinger equation (110j) . it is not difficult 
to obtain the corresponding equations for the new unknown matrices F( m /° ut '(E, r). To 
this end, the ansatz (|2ip is simply substituted into Eq. (fTUl) . After substitution, the 
equation is transformed and simplified, using the following: firstly, the fact that matrices 
WW "*) (E, r) solve Eq. (fl5]l (i.e. Eq. ([10]) without the right-hand side); secondly, the 
Lagrange condition ([22 j) ; thirdly, introducing the diagonal matrix of the channel momenta 



K 



(k x • • 
k 2 ■■ 

•• 







k 



(23) 



nJ 
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and finally using known Wronskian of the Riccati-Hankel functions, 

W {[a) [d r W {out) ] - [d r W {la) ] W {out) = 2iK . (24) 
The derivation can be found in Ref. [2T] . 

As a result, we obtain the following system of first-order differential matrix-equations 
9 r J-(») = -—K- l W {out) V [W^ J" (in) + ^ (out) 7" (out) ] , (25) 

9 r J-(° ut ) = ^k^W^V [W^F^ + W (out) J r{out) ] , (26) 

which are equivalent to the initial second-order Schrodinger equation ( ITOj) . Since these 
equations are of the first order, the boundary conditions for them can only be imposed 
at a single point. A natural way of doing this is to demand that matrix f[2"Tj) is regular at 
the origin. Then 

J- (in) (£,0) = T^iE^O) (27) 

because both (z) and h!f\z) are singular at z = 0, but their singularities exactly 
cancel each other in the combination 

hf\z) + h { -\z) = 2 3f Xz) . (28) 

The choice of the common value for the functions in Eq. (1271) determines the overall 
normalization of the fundamental matrix fl2"Tj) . because Eqs. (I25|I26I) are linear and homo- 
geneous. In the single- channel case, the regular solution is usually normalized in such a 
way that it coincides with the Riccati-Bessel function when r — Y [22]. From Eq. (12"B"|) . it 
is clear that if we demand that our approach gives the traditional single-channel solution 
when iV = 1, we should use the following boundary conditions for Eqs. (I25|26p 

F (in) (E, 0) = ^ out) (E, 0) = - / , (29) 
where I is the diagonal unit matrix. 



Thanks to the fact that the in-coming and out-going waves in the fundamental matrix 
f T2~T]) are factorized, it is easy to construct the physical solutions with proper asymptotics. 
For example, the bound states at large distances can have only the out-going waves. This 
is achieved by finding such combination coefficients in Eq. f|T4|) . i.e. in 



u. 



V \w {in) F (in) 

/ , y V nn"- r n"n' 



, w( out ) -r-( out ) 
VV nn" • r n"n' 



(30) 
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that all the columns involving the in-coming waves, disappear when r — > oo. In other 
words, for the physical wave function of a bound state, we have 




£ F n?(£)o, = o. 



(31) 



n' 



n' 



This is a system of linear homogeneous equations for the unknown combination coefficients 
C n . It has a non-trivial solution if and only if the corresponding determinant is zero. This 
results in Eq. (I2"0~j) which can have solutions at discrete values of the energy. Actually, 
the resonant states are found in the same way, i.e. by solving the same Eq. (I20p . but at 
complex values of S n . This approach therefore offers a unified way of finding both the 
bound and resonant states as well as the S-matrix ffT9l) for the scattering states. The 
corresponding physical wave functions are obtained with correct asymptotic behaviour 
(analytically factorized Riccati-Hankel functions). 

5 Riemann surface 

For a fixed (generally speaking, complex) value of the energy E, each of the N channel 
momenta (j!2p can have two different values 



depending on the choice of the sign in front of the square root. All these momenta are 
involved as papameters in Eqs. ( 1251 1261) . This means that the Jost matrices F( in / out ) (E) 
are not single-valued functions of E. At each point E, they have 2 N different values, for 
all possible combinations of the signs of N channel momenta. 

In complex analysis, the multi-valued functions are treated as single-valued, but de- 
fined on a multi-layered complex surface which is called Riemann surface. In our case, 
each layer (sheet) of this surface corresponds to a different combination of the signs of N 
channel momenta, and thus the Riemann surface of the energy consists of 2 N sheets. 

When we move around a threshold point, we go from one layer to another. Indeed, a 
point on a circle centered at the threshold E n , can be parametrized as E = E n + pexp(itp), 
where p is the distance from E n and (p is the polar angle. The corresponding channel 
momentum 





(32) 




(33) 
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changes its sign after one full circle and comes back to its initial value after two full circles 
(ip = Att). This means that the sheets of the Riemann surface are connected to each other 
and thus form a united multi- layer manifold. The threshold points E n (n = 1, 2, . . . , N) 
are the branching points on this manifold. By moving around these points, we can con- 
tinuously reach any of the 2 N sheets. 

In principle, we can construct the Riemann surface rather arbitrarily by making cuts 
and appropriate connections of the layers. In quantum theory, it is standard that each 
layer is cut along the real energy axis. The cut starts at the branching point and goes 
to infinity in the positive direction. The edges of the cuts of different layers are intercon- 
nected in such a way that the corresponding channel momenta appropriately change their 
signs. 

The simplest two-layer Riemann surface for the single-channel case is easy to visualize 
(see, Fig. [1]). The two-channel problem with two branching points and four intercon- 
nected layers is much more involved. These connections for the three intervals, E < E±, 
E 1 < E < E 2 , and E > E 2 are schematically shown in Fig. [2] [28] - When N > 2, the 
surface is so complicated that it is not worthwhile even trying to visualize it. 

In the present paper, we construct the Jost matrices in such a way that in their 
matrix elements the dependences on odd powers of all channel momenta are factorized 
analytically (see Sec. [7]). The remaining matrices depend only on even powers of all the 
momenta k n and thus are single-valued functions of variable E. This saves us the trouble 
of dealing with the complicated multi-layered manifold. Moreover, using the analytically 
factorized dependence on k n , we can establish some of the symmetry properties of Jost 
matrices, i.e. we can relate their values at some points belonging to different layers of the 
Riemann surface (see Sec. E]). 

6 Complex rotation 

When the potential is cut off at certain radius R, the right-hand sides of Eqs. (I25|26p 
vanish for r > R and the derivatives dj-F^ / ^ become zero, i.e. these functions do not 
change beyond this point. Therefore, in the spirit of the variable phase approach, the 
functions jr(W°ut) ^ 

r) are the Jost matrices for the potential which is cut off at the 
point r. In general, when the potential asymptotically vanishes at large distances, we 
have 




(34) 
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Therefore, the Jost matrices can be calculated by numerical integration of the differential 
equations (I25f26p from r = up to a sufficiently large radius R where the limit (1511) is 
reached within a required accuracy. 

This works perfectly for real values of the energy E. However, when we consider 
complex energies (for example, in search for resonances), a technical difficulty arises. 
This difficulty is caused by the asymptotic behaviour of the Riccati-Hankel functions [2H] , 

h^\kr) — > ^fiexp(±ikr^fi — J . (35) 

As is seen, when k is complex, either hg(kr) or h\ \kr) exponentially diverges, depend- 
ing on the sign of Im k. As a result, either the first or the second of the equations (I25|26p 
does not give a numerically convergent solution. This difficulty is circumvented by using 
the deformed integration path shown in Fig. |3J Instead of integrating the differential 
equations along the real axis from r = to r = R, we can reach the final point via the 
intermediate point r = R' in the complex plane. Moreover, we can safely ignore the arc 
R'R since the potential is practically zero at that distance. 

Why does this complex rotation help? The answer can be found by looking at Eq. 
f )35|) . Indeed, the asymptotic behaviour (divergent or convergent) of the functions h^\kr) 
is determined by the sign of Im(fcr). If k = \k\e lip , we can always find such a rotation 
angle 9 in r = \r\e td that the product 

kr = \kr\e i(v+9) 

has either positive or negative (or even zero) imaginary part. Various technical details of 
using complex rotation in calculating the Jost functions and Jost matrices can be found 
inRefs. [^l EH [321 [331 EH ESI [361 EZl EH [39]. 



7 Factorization 

Eqs. (I25|26p are very convenient for numerical calculation of the Jost matrices. However, 
for the purpose of power-series expansion of these matrices, we need to further transform 
them. The idea of such a transformation is based on the relation between the two pairs 
of linearly independent solutions of the Riccati-Bessel equation, namely, between the 
Riccati-Hankel functions and the pair of Riccati-Bessel je and Riccati-Neumann ye 
functions, 

hf ) {z)=3i(z)±iy t {z) . (36) 
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Introducing the diagonal matrices 



J = - [w {[n) + W {out) ] 



fjhihr) 
ji 2 {k 2 r) 



Y = ~ [W (in) - W {out) ] 



V o 








Vt*(k 2 r) 





ji N (k N r)J 

o \ 



ye N (k N r)J 



(37) 



(3f 



\ o o 

as well as the new unknown matrices 

A(E, r) = J^ (in) (E, r) + J^ (out) (E, r) , (39) 

B(E,r) = i[T {in \E,r)-^ out) (E,r)] , (40) 

we obtain another (equivalent) representation of the fundamental matrix of regular solu- 
tions, 

<$>{E, r) = J(E, r)A(E, r) - Y(E, r)B(E, r) . (41) 
Combining Eqs. (I25|26l) . it is easy to obtain an equivalent system of differential equations 



-K~ 1 YV (J A — YB) 
-K~ l JV (JA-YB) 



for the new unknown matrices, 

d r A = 
d r B = 

with the boundary conditions 

A(E,0) = I, B(E,0) = 

which immediately follow from (|2"9"|) . Similarly to the limit 
converge to their asymptotic values 



(42) 
(43) 

(44) 

these matrices also should 



A{E,r) 



A(E), B(E,r) 



B(E) 



from which the Jost matrices can be obtained, 



E^ in \E) = ~ [A{E) - iB{E)\ , F^\E) = ~ [A(E) + iB{E)\ . 



(45) 



(46) 
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Now, we use the fact that the Riccati-Bessel and Riccati-Neumann functions can be 
represented by absolutely convergent series, 

jl(kr) = ( k z) ,+1 f = r) (47) 



2 J 



Ve(kr) 



n=0 

£ oo 



T(£ + 3/2 + n)n\ \2 J 



kr I ^ TV- 



-1) 



n+i+l 



n=0 



T(-i + 1/2 + n)n\ \ 2 



kr 



2/i 



- =k-%(E,r) 



(46 



where we factorize the the functions and y^, which do not depend on odd powers of k 
and thus are single-valued functions of the energy E. 



Let us look for the matrices A and B in the form 



k^ +1 - 
A.. = -1 A.. 



Bij — k^kj Bij , 



(49) 



where certain powers of the channel momenta are factorized in each individual matrix 
element. When the representations (14"7| |4"51 USD are substituted into Eqs. (I42f430 . all 
the channel- momenta factors cancel out, and we remain with the equations for the tilded 
functions, 



d r A 



d r B 



-Yv(jA-Y&) , 



-JY \JA-YB) 



(50) 



(51) 



where the matrices J and Y differ from fl37j) and fl38l) by the diagonal factors 

(k[^ +1 ■•• \ (k^ ••• \ 

k^ +l ■•• 



J 



V o 







J 



Y 



u-ta 
ft, 2 



V o 







k~ iN 



Y . 



(52) 



The main advantage of Eqs. (150115 ip is that they do not involve any coefficients or functions 
depending on odd powers of the channel momenta. This means that their solutions, i.e. 
the matrices A(E,r) and B(E,r), are single- valued functions of the energy. The multi- 
valuedness of the initial matrices A(E, r) and B(E, r) as well as the fact that they are 
defined on a complicated Riemann surface, are determined by the momentum-factors 
separated in Eqs. (}I9"]) . 
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8 Symmetry of the Jost matrices 

As an example of usefulness of the analytic structure of the Jost matrices, established in 
the previous section, let us consider the relation between and F^ ^. If A(E) and 
B(E) are the asymptotic values of A(E, r) and B(E,r), respectively, then according to 
Eqs. ( H6l H9l) . we have the following semi-analytic expressions for the Jost matrices 

+1 +i 

F (in) ^n a R /roN 

1 mn ~ 2fc£rr+1 mn '~> L 'um • \ ■>■->! 



Ill 



r(out) _ h n a , ^m%_D (ca\ 

n ~ 2k im+1 n 2 " ' 



/it 



If we change the signs of all the channel momenta to the opposite, the matrices A(E) 
and B(E) remain unchanged while the factorized momenta generate a common factor 
(— l) £m+£n = (— \Y m ~ ln and the sign between the two terms in Eqs. ( 1531 E5D is also 
changed to the opposite. In other words, 

FW(-k h -k 2 , . . . , -k N ) = (-lY^F^ih, k 2 ,..., k N ) , (55) 

where for the sake of clarity, we replaced the single independent variable E with the set 
of channel momenta. This means that the two Jost matrices _F < - m - ) and F^ ^ are not 
completely independent. Their values at different points on the Riemann surface are the 
same. The S matrix can be re- written as 

Smn = (-lY^F^i-h, -k 2 , -k N ) [F^ih, k 2 ,..., k N )] _1 , (56) 

which is well known in the case of a single-channel problem. The symmetry property ( 155]) 
is the simplest one (but perhaps the most important) and can be established in a different 
way (see, for example, Ref. [25]). Its simple proof is given here to merely demonstrate 
how the factorized semi-analytic expressions ( 1B^| can be used. Establishing the other 
symmetry properties is beyond the scope of the present paper. 



9 Power-series expansion 

The coefficients of differential equations (I50|5ip . i.e. the products of J(E, r), Y(E, r), and 
the potential matrix V, depend on the parameter E, but do not involve any functions 
depending on channel momenta. The boundary conditions for them, 

A(E,0) = I, B(E,0) = 0, (57) 
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are also momentum independent. This means that the matrices A(E,r) and B(E,r) are 
single-valued functions of the parameter E, i.e. they are defined on a single sheet of the 
complex E-p\a.ne. 

Since the matrices J(E, r) and Y(E, r) are holomorphic, the coefficients of differential 
equations (I50|I51I) are also holomorphic for any finite radius r. Moreover, the boundary 
conditions (|57() do not depend on E. According to the Poincare theorem [3D], this means 
that A(E, r) and B(E, r) are holomorphic functions of the parameter E for any finite 
value of r. Therefore they can be expanded in the power series around an arbitrary point 
Eq on the complex plane of the energy, 

oo 

A(E,r) = ^(£-£ ) n a„(£o,r) , (58) 

n=0 

oo 

B(E,r) = Y,( E - E o) n Pn(E ,r) , (59) 

n=0 

where the unknown expansion coefficients a n and /3 n are (N x iV)-matrices depending not 
only on variable r but also on the choice of the point Eq. 

Eqs. (I50|l5ip involve matrices J and Y for which we can obtain the expansions of the 
same kind but with known coefficients, 7 n and r) n , 

oo 

J(E,r) = J2( E - E o) n ln(E ,r) , (60) 

n=0 

oo 

Y(E,r) = ^(E-E o r Vn (E ,r) . (61) 

n=0 

Simple recurrency relations for calculating the expansion coefficients (matrices) 7 n (£'o,r) 
and T] n (E ,r) are derived in the Appendix IA1 

Substituting the expansions fl58| |59| [60] IBTj) into Eqs. f l5"0f5~T|) . and equating the 
coefficients of the same powers of (E — Eq), we obtain the following system of differential 
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equations for the unknown matrices a n and /3 n 

d r a n - - ^ ViV{lj®k ~ VjPk) , (62) 

i+j+k=n 

d r (3n = - 7^ 7i^(7j"fc - , n = 0, 1, 2, . . . (63) 

The boundary conditions ( 15"T|) are independent of E. As is easy to see, this implies that 

a n (E Q ,0) = 5 n0 I , /3 n (E ,0) = 0. (64) 

In other words, all these matrices vanish at the origin, except for the matrix a which 
becomes a diagonal unit matrix at r = 0. 

Therefore for any finite radius r and an arbitrary complex Eq, we have a simple proce- 
dure for calculating the expansion coefficients. What we actually need are the correspond- 
ing expansions of the matrices A(E,r) and B(E,r) when r — > oo, from which the Jost 
matrices fl4"6"j) can be constructed. In this context, taking the limit r — > oo is not a trivial 
procedure. Indeed, at large distances the Riccati-Bessel and Riccati-Neumann functions 
behave as linear combinations of the exponential functions (|35|) . There is nothing wrong 
in this behaviour if the channel momenta are real. If however the energy is complex and 
the momenta have non-zero imaginary parts, then the matrices J(E,r) and Y(E,r) in- 
volve divergent exponential functions and thus are not holomorphic when r — > oo. 

To some extent the situation can be saved by using exponentially decaying poten- 
tials V nn i{r) ~ exp(— A nn /r), which compensate the divergence of J(E,r) and Y(E,r) 
within certain domain T> of the complex .E-plane along its real axis. The borders of 
the domain T> are determined by the requirement that none of the coefficients of Eqs. 
f )5T?|5T]) are divergent. The behaviour of these coefficients is determined by the prod- 
ucts exp(±ik n r) exp(— A nn /r) exp(±i£vr). The worst case is when we have the product 
of the most rapidly growing exponentials stemming from the Riccati-Bessel and Riccati- 
Neumann functions, and the most slowly decaying element of the potential matrix. There- 
fore 



V 



\ E : max 21m ^Jl[i n {E — E n )/h 2 < min \ nn > \ . (65) 



The faster the potential matrix decays, the wider is the domain. An example of such a 
domain is shown in Fig. H]for the two-channel model used in Sec. [10J 
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Therefore when r — > oo, we can use the Poincare theorem only within D. This means 
that only within this domain the matrices A(E, oo) and B(E, oo) are holomorphic func- 
tions of the parameter E. A rigorous analysis (using a different method) of the analyticity 
domain of the multi-channel T-matrix for the same class of potentials is given by Mo- 
tovilov in Ref. |23|. 



If the potential matrix V(r) (or at least its long-range tail) is an analytic function of 
complex variable r and exponentially decays along any ray r = \r\e ld within certain sector 
#min < 6 < # max of the complex r-plane, then the domain of analyticity of the matrices 
A(E, oo) and B(E, oo) can be extended by using the complex rotation described in Sec. 

The physically interesting domain of the E-pl&ne where the expansion proposed in the 
present paper, can be used in practical calculations, lies on the positive real axis (scatter- 
ing and reactions) and in the close vicinity below it (pronounced resonances). Therefore 
we can say that to all practical purposes the expansion can be done near an arbitrary 
point E . 



If we denote the asymptotic values of the matrices a n and (3 n as 

a n (E ,r) — > a n (E ) , and (3 n (E ,r) — > b n (E ) 



(66) 



then in vicinity of any chosen point E on the Riemann surface, we can obtain semi- 
analytic expressions for the Jost matrices in the form 



M 



3=0 



K n 



2k e ™ +1 y<lj > mn 



j /mn 



(67) 



M 



p(out) 



3=0 



E n V 



h.tn+1 

2^ +1 



j /mn 



(68) 



where M is the maximal power of our expansion. The non-analytic quantities in Eqs. 
<E7\ IB"8"|) . i.e. the matrices a,j and bj, are obtained by numerical integration of differential 



equations ( 16"2"1 [B^]) with the boundary conditions from r = to certain large value 
r = R (in some cases, this should be done along the deformed contour shown in Fig. [3]). 
For a given energy E, the choice of the sheet of the Riemann surface, where the Jost 
matrices are considered, is done by appropriately choosing the signs in front of the square 
roots fl32l) for calculating the channel momenta used in Eqs. ( 16?] |68|) . The central point 
i?o of the expansion as well as the numerically obtained matrices dj and bj are the same 
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for all the layers of the Riemann surface. 

In principle, the power-series expansions fl67] 168]) includes infinite number of terms. 
In practice, however, we may take into account just a few terms (M < oo) and therefore 
have to solve certain number of equations in the infinite system fl62l [631) . It should be 
emphasized that not all equations of this system are coupled to each other. Indeed, due 
to the condition i + j + k = n on their right-hand sides, the equations for % and (3m (for 
any M ^ 0) are only linked to the corresponding equations with n < M. For example, 
the first pair of equations, 

d r a = -r] V(-f a - r] (3 ) , (69) 

d r /3o = -70^(70^0 - VaPo) , (70) 

is self-contained and is not linked to any other equation of the system. The second pair 
is linked only to the first one, and so on. This means that by considering a finite number 
of equations of this system, we do not introduce a truncation error. 

In principle, it is possible to find some physical meaning (scattering length, effective 
radius etc.) of certain combinations of the expansion coefficients (cij) mn and (bj) mn in 
Eqs. f)67|68l) . This would require the parametrization of the corresponding S- matrix f fT9|) 
in terms of the channel phase-shifts. This is easily done for a single channel problem (see 
Ref. [16J). However for a multi-channel system (even with minimal N = 2) this results 
in rather complicated formulae which do not add any clarity to the method. Moreover, 
the main advantage of the proposed method consists in the fact that it can be used not 
only at low but at any energies (even complex) where the notions of scattering length and 
effective radius are meaningless anyway. This is why we think that within our method 
it is not worthwhile to derive any expressions for something like multi-channel scattering 
length and effective radius. 



10 Numerical example 

In order to demonstrate the efficiency and accuracy of the proposed method, we do nu- 
merical calculations for a well-studied model. For this purpose, we use the two-channel 
potential suggested by Noro and Taylor 




V(r) = r 2 e~ r , (71) 
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which is given in arbitrary units such that fix = H2 = he = 1. The threshold energies for 
the two channels of the Noro and Taylor model are E\ = and E 2 = 0.1, and the angular 
momentum is zero in both channels, l\ = £2 = 0. This potential has an attractive well 
in the first channel, a repulsive barrier in the second, and rather strong coupling between 
the channels. As a result, it generates a rich spectrum of bound and resonant states (see 
Fig. \5§ as well as a non-trivial energy dependence of the channel and transition cross 
sections (see Figs. El El and[S]) [27J. This model is therefore a difficult testing ground for 
any new method designed for describing multi-channel processes. 

As a first test, we do the power series expansions of the Jost matrices ( 1671 [68]) around 
the point E = 5 + iO on the real axis and with six terms, i.e. with M = 5. The central 
point of the expansion was chosen to be not far from the first resonance, where the channel 
cross sections have some non-trivial energy dependence. As is seen in Figs. [61 [7J and El 
the thin curves representing the approximate cross sections, reproduce the corresponding 
exact cross sections rather well within a wide energy interval. 

Since the approximate Jost matrices flBTl |6"8|) coincide with the exact matrices on a 
segment of the real axis, they must be valid also at the nearby points of the complex 
energy surface. Comparing the exact and approximate values of det (E) at complex 

E around the point E on the third sheet ( ) of the Riemann surface (see Fig. [2]), we 

found the domains within which the relative accuracy is better than 1%, 5%, and 10%. 
These domains are shown in Fig. 

It is seen that the first resonance is within the domain of 1% accuracy and there- 
fore must be reproduced by the approximate Jost matrix. And indeed, its determi- 
nant has zero at E — 4.768178 — iO. 000686 which is very close to the exact value 
E = 4.768197 — iO. 000710. With more terms in the expansion, the difference becomes 
smaller, and with M — 10 all the digits are the same. 

The other place where we tested the expansion, was the point Eq — 7.5 — z2.0 on the 

third sheet ( ) of the Riemann surface. This point is almost in the middle between 

the second and the third resonances. Fig. [TU] shows how the domain of 1% accuracy 
increases with increasing number of terms (M = 5, 7, 10) in the expansions (I6T1 I6"81) . As 
is seen in Fig. [101 even with M = 5 both resonances are reproduced relatively well (the 
filled and open circles represent the exact and approximate positions of the resonances). 
The zeros of the exact Jost matrix determinat are at E — 7.241200 — iO. 755956 and 
E = 8.171217-i3.254166 while the expansion with M = 5 gives E = 7.131204-^0.768670 
and E = 8.241795 — i2. 982867. We deliberately chose the point E far away from both 



20 



resonances. If it is close to any of them, the resonance can be found very accurately 

11 Conclusion 

In the present paper, we show that each matrix element of multi-channel Jost matrix can 
be written as a sum of two terms, and each term can be factorized in such a way that it 
assumes the form of a product of certain combination of the channel momenta k n times an 
analytic single-valued function of the energy E. This means that all the branching points 
of the Riemann energy-surface are given in the Jost matrix explicitly via the channel- 
momentum factors. The remaining energy-dependent factors in all its matrix elements 
are defined on single energy plane which does not have any branching points anymore. 

For these energy-dependent functions, we derive a system of first-order differential 
equations. Then, using the fact that the functions are analytic, we expand them in the 
power series and obtain a system of differential equations that determine the expansion 
coefficients. A systematic procedure developed in the present paper, allows us to accu- 
rately calculate the power-series expansion of the Jost matrices practically at any point on 
the Riemann surface of the energy. Actually, the expansion is done for the single-valued 
functions of the energy, while the choice of the sheet of the Riemann surface is done by 
appropriately choosing the signs of the channel momenta in the momentum-dependent 
factors. 

The expansion suggested in the present paper, makes it possible to obtain a semi- 
analytic expression for the Jost-matrix (and therefore for the S-matrix) near an arbitrary 
point on the Riemann surface and thus to locate the spectral points (bound and resonant 
states) as the S-matrix poles. Alternatively, the expansion can be used to parametrize 
experimental data, where the unknown expansion coefficients are the fitting parameters. 
Such a parametrization will have the correct analytic properties. After fitting the data 
given at real energies, one can use the semi-analytic Jost matrix to search for resonances in 
the nearby domain of the Riemann surface. The efficiency and accuracy of the suggested 
expansion is demonstrated by an example of a two-channel model. 

APPENDIX: 
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A Expansion of the Riccati functions 

As is given by Eqs. (1471 |4"8"|) . the Riccati-Bessel and Riccati- Neumann functions can be 
factorized as 

j e (kr) = k e+1 ] e (E, r) , y e (kr) = k~%(E, r) , (72) 

where the tilded functions depend on k 2 , i.e. on the energy. These functions are holomor- 
phic and thus can be expanded in Taylor series at an arbitrary point E = E , 



Jt 



i £ {E, r) = J2(E - E ) n g in (E , r) , y e (E, r) = ^(E - E ) n t in {E , 



n=0 



n=0 



where the coefficients are given by the derivatives 



g en (E ,r) 



nl 



h(E,r) 



t en (E ,r) = — 



dE 
d" 



dE r 



Ve(E,r) 



E=E 



E=E 



1 

n\ 
1 

nl 



d n je(kr) 
dE n k e+1 

d" 



dE' 



: k ye(kr) 



E=E 



E=E 



In order to find these derivatives, we use the following relations [29J 



d 

dz 



Jew 



Je+i{z) 
i+i 



and 



z 



d 

dz 



[z e y E (z)] = /yi-i{z) . 



After a simple but lengthy algebra, we finally obtain 

g in (E ,r) 



1 




| n ]i+n{kr)~ 




n\ 




1 fcl+n+l 


E=E 



ten(E ,r) = ^(^) n [k^ n (kr)} E ^ Eo 



(73) 



(74) 
(75) 

(76) 



(77) 
(78) 



The matrices 7 n and r) n of Eqs. fl60l IBTj) are diagonal with each row having the functions 
( 1771 1751) with ji and k for the corresponding channel. These functions should be the same 
for all sheets of the Riemann surface, i.e. for any choice of the signs of channel momenta. 
This is so indeed since jt(—z) = (— lY +1 ji(z) and yi(—z) = (—l) e yi(z). 
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68.503722 
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Table 1: Spectral points E = E r — iT/2 generated by the potential ( 17T1) and shown in 
Fig. El 
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Figure 1: Riemann surface of the energy for a single-channel problem. 



(++) 

(-+) 

( — ) 

(+-) 

E < E x 

Figure 2: Schematically shown interconnections of the layers of the Riemann surface for a 
two-channel problem at three different energy intervals. The layers correspond to different 
combinations of the signs (indicated in brackets) of Im/ci and ImA^. 
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Figure 3: A deformed path for integrating the differential equations. 
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Figure 4: The domain V defined by Eq. ( 1551) , for the two-channel model described in Sec. 
[TU1 The matrices A(E, oo) and B(E, oo) for this model are holomorphic in the area enclosed 
between the parabolic curves. The parabola crosses the real axis at E = —0.025. 
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Figure 5: Spectral points generated by the potential (TTTT) and given in Table [U 
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S-matrix poles 
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Figure 6: Energy dependence of the elastic scattering cross section in channel 1 for the 
potential (17TI) . Few of the ^-matrix poles (see Tableland Fig. [5]) are shown in the lower part 
of the Figure. The thick curve represents the exact cross section, while the thin curve shows 
the cross section obtained with the expansions (1571 loll where E = 5 + iO and M = 5. 



30 



<r(l 2) 
0.5 - 
0.4 - 
0.3 - 
0.2 - 



expansion 



0.05 




i i r i i i i i i | 
4.765 4.775 




Figure 7: Cross section energy dependence of the inelastic transition (1 — > 2) for the potential 
(I7T1) . The thick curve represents the exact cross section, while the thin curve shows the cross 
section obtained with the expansions (!o7l 1551) where E = 5 + iO and M = 5. 
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Figure 8: Energy dependence of the elastic scattering cross section in channel 2 for the 
potential (I7TI) . The thick curve represents the exact cross section, while the thin curve shows 
the cross section obtained with the expansions (loTl 1551) where E = 5 + iO and M = 5. 
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Figure 9: The domains within which the Jost matrix determinant for the potential ( 1711 ) is 
reproduced, using the first five terms (M = 5) of the expansion (!67f) , with the accuracy better 
than 1%, 5% and 10%. The expansion was done around the point E = 5 on the real axis. 
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Figure 10: Growth of the 1% accuracy domain with the increase of the number M of terms in 
expansions ( !67i ). which were done around E = 7.5 — i2.0 for the potential ( I7T1 ). Filled circles 
indicate the exact position of two resonances, while the open circles are their approximate 
positions obtained with M — 5. 
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